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ABSTRACT 

Early-type dwarfs are the most common galaxy in the local universe, yet their origin and 
evolution remain a mystery. Various cosmological scenarios predict that dwarf-like galaxies in 
dense areas are the first to form and hence should be the oldest stellar systems in clusters. 
By using radial velocities of early-type dwarfs in the Virgo cluster we demonstrate that these 
galaxies are not an old cluster population but have signatures of production from the infall of 
field galaxies. Evidence of this includes the combined large dispersions and sub-structure in 
spatial and kinematic distributions for Virgo early-type dwarfs, and a velocity dispersion ratio 
with giant ellipticals expected for virialized and accreted populations. We also argue that these 
galaxies cannot originate from accreted field dwarfs, but must have physically evolved from a 
precursor population, of different morphology, that fell into Virgo some time in the past. 



Subject headings: Galaxy Clusters: Virgo, Galaxies: Dwarf Ellipticals/Spheroidals, Galaxies: 
Formation, Evolution 



1. Introduction 

Early-type dwarfs are the most common galaxies in the local universe, yet we know very little about 
them (cf. Gallagher & Wyse 1994; Ferguson & Binggeli 1994). This is due to their very low luminosities 
and surface brightnesses hindering any detailed studies. Luminosity functions of nearby clusters show 
that faint galaxies^ are by far more common than brighter galaxies (e.g., de Propris et al. 1995; Seeker, 
Harris & Plummer 1997; Trentham 1997; Phillipps et al. 1998). Additionally, in comparison to the field, 
galaxy clusters have steeper luminosity functions dominated by dE-like galaxies (e.g., Ferguson & Sandage 
1989; de Propris et al. 1995 Seeker et al. 1997; Marzke & da Costa 1997; Phillipps et al. 1998; Mobasher 
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& Trentham 1998). How these objects originated, and continue to exist in large numbers in clusters is 
presently a mystery. 

A very basic question to ask is: where did these cluster dwarf galaxies come from? Do they have 
an origin similar to the dwarf elliptical and spheroidal galaxies seen in groups such as the Local Group 
(cf. Mateo 1998; van den Bergh 2000)? One basic difference between group and cluster dwarfs is the 
way they are distributed. In groups of galaxies, dwarf spheroidals are for the most part clustered around 
giants (Mateo 1998; Grebel 1999). However, in clusters dwarfs are in general not distributed around the 
giants, but follow the general potential (e.g., Binggeli, Tarenghi & Sandage 1990; Ferguson 1992; Seeker 
et al. 1997). In other respects, such as surface brightness, magnitude, and color, bright dEs in clusters 
have properties similar to bright Local Group dEs (Bothun et al. 1985; James 1991; Held & Mould 1994; 
Phillipps et al. 1998). 

Dwarfs in rich clusters are faint and hard to identify and as a result they were not discovered until 
Reeves (1956, 1977) and Hodge (1959, 1960, 1965) studied them in the Virgo and Fornax clusters. Since 
then many observations, conducted at various wavelengths and employing various techniques, have been 
applied to cluster dwarfs, with particular activity in the 1980s. Many of these studies were mainly interested 
in determining the galaxy luminosity functions in clusters (Ferguson & Sandage 1988; Phillipps et al. 1998; 
Mobasher & Trentham 1998; Trentham 1998) and not necessarily in the properties of the dwarf galaxies 
themselves. Luminosity functions are important for cosmology and galaxy formation theories, but they also 
reveal what the relative number densities of cluster dwarfs are in comparison with the field. 

There arc at least four different ideas for explaining where dEs seen in Virgo and other galaxy clusters 
originate. These are: (1) dEs are old primordial galaxies that formed early in the Virgo proto-cluster (e.g., 
White & Frenk 1991; Ostriker 1993; Farooqui 1994) and have survived in the cluster environment for several 
Gyrs. (2) dEs are the passively evolved compact blue and emission line galaxies seen at redshifts z <~ 1 
(e.g., Babul & Rees 1992; Koo et al. 1995; Koo et al. 1997) (3) dEs are stripped and 'harassed' spirals 
accreted into clusters, creating a morphological transformation: spiral — > dE (Mao & Mo 1998; Moore et 
al. 1998; Moore et al. 1999) (4) dEs evolved from dwarf irregulars (Lin & Faber 1983). Alternatively dwarf 
ellipticals may form though a combination of these scenarios. 

Solution (4) was considered quite seriously and given much attention in the 1980s, but has since been 
largely discarded (e.g., Bothun et al. 1986; Gallagher & Hunter 1989) due to the problem of explaining 
large populations of dEs that are more massive, brighter, and more compact than typical dwarf irregulars. 
Since any evolution from dlrr to dE requires fading, bright dEs cannot be accounted for by this evolution. 

Possibility (2) was inspired by the large numbers of faint blue galaxies seen in rcdshift surveys (Babul 
& Rees 1992). These blue galaxies have low masses, similar to masses of dwarfs (Koo et al. 1995). If 
no further star formation occurs then the magnitudes and surface brightnesses of these objects would be 
similar to those of dwarf ellipticals (Koo et al. 1995), thus allowing for the possibility of an evolutionary 
connection. However, there are several potential problems with this idea, the most important being the fact 
that many Local Group dwarf ellipticals contain several episodes of star formation (cf. Mateo 1998) and 
hence have not evolved passively. A second problem is that the faint blue galaxies are a field population, 
yet in the nearby universe most dEs are in clusters. Third, some Virgo dEs are brighter than the luminosity 
to which these high redshift, star-forming galaxies would evolve with no further star formation. Thus, a 
strictly passive evolution from a burst of star formation at z <~ 1 to the present cannot account for cluster 
dEs. 

In this paper, we mainly address, and try to distinguish between, possibilities (1) and (3). The 
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hypothesis we test is that cluster dEs form either by the galaxy harassment scenario (e.g., Moore et al. 
1998), or are an old cluster population. These models depend upon a very active source of accreted galaxies 
to account for the large population of dEs in clusters. For the harassment idea to be correct, there must 
be signatures that dEs were accreted and underwent mass stripping, or that a large population of spirals 
existed in the proto-cluster. Modern Virgo spirals have properties that suggest they have recently been 
accreted (Huchra 1984) and this process likely occurred in the past as well. Harassment has also been used 
to explain the morphological transformation of Butcher-Oemler galaxies seen at moderate cluster redshifts 
(z«0.8) into modern dwarf galaxies (Moore ct al. 1996). A galaxy in a cluster, such as Virgo, interacts 
weakly with the other cluster members due to the high relative velocities between systems. Stars in the 
outer parts of accreted galaxies are stripped away due to the increase of internal energy from impulsive 
interactions and at the same time gas is stripped by ram pressure (Lee, Richer, & McCall 2000; Mori & 
Burkert 2000). Over the 1 Gyr orbital time the galaxy begins to lose its stars and gas, and over a longer 
time can become a spheroidal system (Moore ct al. 1998). Before a galaxy is stripped any gas it contains 
can become compressed by cluster tides producing a starburst. Harassment does not differentiate between 
accreted galaxies and cluster members or galaxies on isotropic or anisotropic orbits. This theory predicts 
that the larger spiral galaxies become SOs (Moore et al. 1999) while the smaller ones will eventually 
resemble dwarf ellipticals. It remains to be determined, but seems unlikely, whether it is possible that each 
dE originated from a spiral. The increase in number of star forming galaxies in z^ 1 clusters may not need 
an explanation from an increase in the infall rate, since the field is rich in starbursts at this same redshift 
(e.g., Sullivan et al. 2000). 

One way to test if cluster dEs originated from a relatively recently accreted population is to determine 
how their velocity and spatial distributions compare to those of giant ellipticals and spirals. The one 
dimensional velocity dispersions of clusters are typically less than 1000 km s _1 (for Virgo the value is 
about 700 km s _1 ). From radial velocity studies of Virgo galaxies it is known that the dispersion of the 
spiral population is higher and has a flatter distribution (less relaxed) than that of the ellipticals which is 
Gaussian, an expected form for a relaxed population (e.g., Huchra 1985, Schindler, Binggeli, & Bohringer 
1999 (SBB99)). This velocity distribution difference is usually interpreted as a consequence of the infall of 
spiral galaxies. As a very basic test we can measure the velocity and spatial distributions of Virgo dwarf 
galaxies to determine if they have such accretion signatures, or instead signatures of an older relaxed cluster 
population. Furthermore, we can compare these distributions to various predictions, including the effects of 
relaxation, virialization, equipartition and dynamical friction. 

We show that the early type dwarfs in Virgo, as well as all Virgo galaxies except ellipticals, have no 
strong signs of relaxation. Based on both their velocity and spatial distributions, we demonstrate that most 
Virgo dEs cannot have existed in Virgo for as long as the giant ellipticals have. The kinematic properties of 
Virgo dwarf ellipticals also suggest that these galaxies originate from an accreted component, whose infall 
rate peaked some time in the past. This paper is organized as follows: §2 describes the observations and 
data, §3 gives our results, §4 describes various dynamical models, §5 is a discussion of our results, and 
finally §6 gives a summary. A Virgo cluster distance of 18 Mpc is used throughout, where 6° = 1.9 Mpc. 

2. Data and Observations 

The Virgo Cluster of galaxies is the nearest rich cluster, presenting both advantages and problems for 
studying its galaxies and structure. It is close enough so that its galaxies can be studied in detail, but 
because it is so close, it is spread out over 10 degrees in the sky, making multiplexed observations of its 
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galaxies time consuming. 

Measuring radial velocities of faint, low surface brightness dwarf galaxies is also technically very 
challenging. The brightest dwarfs can be studied with little difficulty, but the fainter ones, with surface 
brightnesses significantly fainter than the night sky, though very numerous in Virgo, are difficult to observe 
spectroscopically. As such we can only determine radial velocities of relatively bright dEs. Previous to this 
work, the last published radial velocity catalog of Virgo contained 91 dE + dSO galaxies (Binggeli, Popescu 
& Tammann 1993). We increase the total number to 142 by measuring new velocities and combining these 
with previously published values. 

We undertook radial velocity observations of selected Virgo galaxies with the Hydra multi-object 
spectrograph on the WIYN 3.5m telescope, and with the RC long-slit spectrograph on the Kitt Peak 
National Observatory 4m telescope. 

Hydra contained 95 fibers at the time these observations were made in 1999 April and 2000 March. 
For these observations Hydra's fibers fed the bench spectrograph with a 400 lines/mm grating during 1999 
April and the 600 lines/mm grating in 2000 March. Total dispersions of 2.1 A pixel -1 and 1.4 A pixel -1 
and equivalent resolutions of 6.9 A and 4.6 A were obtained for the two gratings respectively, giving a 
typical error of 35 - 50 km s -1 in radial velocity. The wavelength ranges for the April 1999 and March 2000 
runs were AA3400 - 7600 A and AA3900 - 6700 A. In both runs we used the 3" diameter Hydra blue fibers. 
The field of view of Hydra is approximately 1° . 

We observed a total of three fields over the course of three nights in 1999 and the same number of 
fields in 2000. Radial velocities were obtained in three different areas of the cluster; roughly the central 
part (near M87) as well as the northern (5 = 13° : near M86) and southern clumps (5 = 7°: near M49). 
Sample selection was done using digitized Palomar sky survey plates, in conjunction with morphological 
classifications from the Binggeli, Sandage & Tammann (1985) Virgo Cluster Catalog (hereafter VCC). 
Each fiber setup contained between 15 to 25 galaxies, with most of the remaining fibers positioned 
on the sky across each field. Some Hydra fibers which could not be assigned to dEs were placed on 
non-dwarf Virgo galaxies whose redshifts were unknown. Total exposure time for each field was typically 
about 11,000 seconds. We bias subtracted, flat-fielded and then wavelength-calibrated the spectra with a 
Helium-Neon-Argon lamp. All reduction procedures, except the basic CCD processing, occurred using the 
IRAF task DOHYDRA. We then combined each of the individual frames and later sky-subtracted galaxy 
fibers using the closest sky fibers to each object. Sky subtraction was typically successful with an effective 
removal of the sky continuum and telluric emission features, with occasional residual sky lines in the red 
part of the spectrum. These remaining faint sky lines, typically in the region near 7000 A, were removed by 
hand. 

Long slit data from the KPNO 4m telescope were used to determine the radial velocities of two dwarf 
ellipticals in June 1999. The RC spectrograph with the KP-007 grating was used giving a dispersion of 
1.4 A pixel -1 , and a wavelength coverage AA 4400 - 7350 A. These spectra were reduced and extracted using 
the APALL program in the IRAF long-slit package. 

During these observations we obtained spectroscopy of velocity standards to calibrate the velocity zero 
point. The stars observed are listed in Table 1; most are G-K spectral types. Each stellar spectrum was 
reduced and calibrated the same way as were the galaxy spectra. For each night we create a single velocity 
template by combining the spectrum of each star observed, using the task SUMSPEC in the RVSAO 
(Radial Velocity Smithsonian Astrophysical Observatory) package. 
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Using these standards we then computed the radial velocities for our complete sample of observed 
objects. Only a fraction, ~ 25%, of the observed galaxies had high enough signal-to-noise (S/N) ratios for a 
reliable radial velocity measurement. Figure 1 shows a typical Hydra dE spectrum. We determine velocities 
by cross-correlating absorption features with the velocity template using the XCSAO task in the RVSAO 
2.0 package within the IRAF environment (see Kurtz & Douglas 1998). XCSAO outputs a radial velocity 
for each object computed from cross correlation fits, including a quality R-value (Tonry & Davis 1979) 
where a higher value of R indicates a better fit, and we adopt R > 3 as the limit of reliability (Kurtz & 
Douglass 1998). Of those observed, 43 have high enough S/N and R values > 3 for a reliable radial velocity 
measurement, and 37 of these are dwarf ellipticals or dwarf lenticulars. In Figure 2, we plot the difference 
between our computed velocities with published radial velocities for those few galaxies with previously 
measured values. The average absolute difference between our measured velocities and those previously 
published is 128±71 km s _1 . All new velocities are shown in Table 2, while Table 3 lists all Virgo dEs with 
known radial velocities. 

We combine our new velocity measurements with previously published ones. The most recent previous 
studies on the dynamical properties of the Virgo cluster, Binggeli et al. (1993) and SBB99 contained 385 
and 403 velocities respectively. The present study uses 497 velocities, with new velocities mainly of faint 
galaxies, taken from a variety of sources including Drinkwater et al. (1996) and Grogin, Geller & Huchra 
(1998). The distribution across Virgo of all the galaxies used in our analysis is shown in Figure 3. 

2.1. Galaxy Morphology and Populations 

Identifying a morphological type for a galaxy, or for our purposes, identifying members of distinct 
populations, is never a simple or straightforward procedure. In Virgo the problem is even greater due to 
the poor correlation of Hubble types with global properties such as star formation rates (Koopmann & 
Kenney 1998). A morphological classification therefore has its limitations. None-the-less, we characterize 7 
distinct galaxy populations for this study based on morphology alone. These are: classical ellipticals (E), 
lenticulars (SO), spirals (Sa-Sd), irregulars (Irr, I, dlrr, Sm, blue compact dwarfs), dwarf ellipticals (dE, 
dSph), nucleated dwarf ellipticals (dE, N), and dwarf lenticulars (dSO). We often combine dwarf ellipticals 
and dwarf nucleated ellipticals into one class which we denote as the total dwarf ellipticals (dET) . 

Almost all of the morphological types are taken from the Virgo Cluster Catalog (VCC). Some galaxies 
in this paper are unclassified in the VCC catalog, particularly in the region 12 h 50 m 00 ;s - 13 /t 15 m 00 s . 
We either classified these galaxies ourselves, using second generation Palomar Sky Survey images, or 
adopted classifications from other published studies. If there were a significant uncertainty in a galaxy's 
morphological type, we assigned no classification and the galaxy was left out of the analysis involving 
galaxy populations. The morphological identifications of Virgo core galaxies (defined in the next section) 
with radial velocities breaks down as follows: 40 ellipticals, 56 SOs, 18 dEs, 79 nucleated dEs, 28 dSOs, 81 
irregulars and 119 spirals. 



3. Results 
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3.1. Virgo Cluster: Structure 

Smith (1936) was the first to study the dynamical structure of the Virgo cluster based on 32 radial 
velocities obtained by himself, V.M. Slipher and M.L. Humason (unpublished). Using this data, Smith 
(1936) concluded that the Virgo cluster is a bound system, and derived a mass for it using the kinematics of 
these galaxies. This provided some of the first evidence for extensive dark matter in clusters, similar to the 
pioneering study of Coma done by Zwicky (1933). This derived mass assumed that the cluster components 
are virialized; however the Virgo cluster is not as indicated by the presence of sub-clusters which have been 
known since the study of de Vaucouleurs (1961). More detailed observations revealed Virgo to be a complex 
system of different aggregate clouds of galaxies that are part of the Local Supercluster (e.g., de Vaucouleurs 
1961; Holmberg 1961; de Vaucouleurs & de Vaucouleurs 1973; Huchra 1985; Binggeli et al. 1987; Binggeli 
ct al. 1993; SBB99). De Vaucouleurs & de Vaucouleurs (1973) identify and named several of these clouds, 
in the area of the Virgo cluster studied here {a (J2000): 12 h - 13 h , 8 (J2000): 0° - 20° ; v < 3000 km s" 1 ). 
These are, using the de Vaucouleurs & de Vaucouleurs (1973) notation, the M cloud, Wa and Wb clouds 
and the Virgo I and II clusters. The W cloud (combined Wa and Wb) is the most distinct, with a well 
defined locus at 5 <~ 5° and velocities near 2600 km s _1 . The Virgo II cluster (the B cluster in Binggeli et 
al. 1987) is south of the area used in this study, but it and the other galaxy clouds near Virgo are probably 
merging with each other based on their relative positions and velocities. 

These various clouds can bias results; thus to study the Virgo core cloud we define two separate 
regions. We define the "Virgo Region" as those galaxies with a (J2000): 12 h - 13'\ S (J2000): 0°- 
20° and with velocities v < 3000 km s _1 . The distribution of velocities for each of our defined galaxy 
populations is plotted in Figure 4. This is still a rather generous spatial and velocity coverage that is 
certainly unsatisfactory for studying only galaxies associated with the primary Virgo core region. To study 
only galaxies in the cluster proper we identify a separate "Virgo Core" containing only galaxies within 6° of 
the center of the luminosity density of the Virgo I or A main cluster (defined by Binggeli ct al. (1987) as 
(B1950) 12 h 27 m 50 s ; +13°01'24") and with radial velocities v< 2400 km s^ 1 . Spatially these galaxies 
are found within the circle shown in Figure 3. Isolation of this sample allows us to perform statistical 
tests without a large bias, or contamination from Local Supercluster clouds seen in projection towards the 
cluster, although some contamination is inevitable. The "Virgo Core" velocity histograms are plotted in 
Figure 5; while the distribution of spatial positions of all likely members of Virgo, based on VCC estimates, 
are shown for each galaxy population in Figure 6. Figure 6 contains an essentially complete listing of all 
Virgo galaxies down to magnitude B — 18. 

The analysis presented here is only concerned with these two areas in Virgo, but mostly with the Virgo 
Core region. While our large sample of 497 radial velocities, 429 of which are in the Virgo Core, permits a 
detailed dynamical study of the cluster and its neighboring clouds, we postpone an analysis of these clouds 
for another paper. It is sufficient for comparison purposes to state a few basic parameters (Table 4). The 
mean Virgo Core heliocentric velocity is found to be 1064±34 km s _1 , less than 0.5tr from the value of 
1050±35 km s _1 found by Binggeli et al. (1993) and has a total velocity dispersion of 705 km s _1 , close to 
the 699 km s^ 1 dispersion found by Binggeli et al. (1993). 

The usefulness of these mean and dispersion quantities is limited since as discussed above, the Virgo 
cluster is not a simple relaxed system with a Gaussian velocity distribution (see Figures 4-8). Figures 7 
shows the velocity of each Virgo galaxy relative to the mean Virgo velocity of 1064 km s _1 , as a function of 
the absolute magnitude. The dashed lines show the escape velocity from the Virgo core computed from the 
mass profile of Giraud (1999). From Figure 7 we see no equipartition of energy, nor any obvious dynamical 
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friction effects that would produce lower velocities for higher mass objects. 

3.1.1. Virgo Clouds: Galaxies Bound to the Cluster? 

Figure 3 shows the velocities of objects across the Virgo cluster. Triangles represent galaxies with 
velocities less than 1000 km s" 1 , with larger symbols signifying larger velocity deviations from 1000 km s -1 . 
Circles are galaxies with velocities larger than 1000 km s _1 , with the largest circles signifying galaxies at v 
> 2000 km s _1 . Signatures of various 'sub-clusters' or clouds can be seen from the different symbols. The 
spatial clustering in Figure 3 is however partially the result of selection biases introduced by spectroscopic 
sample selection. 

There is evidence for spatial clustering in Virgo, as is shown in Figure 6. SBB99 found that galaxy 
clustering correlates with the intensity of X-ray emitting gas in Virgo. In the north, there is a sub-cluster 
with lower velocities associated with the Virgo cloud around M86. Likewise, we see velocity structures near 
M87 (the X in Figure 3) and the W cloud near (J2000) (a, 6) w (12.35, 8.00). 

The dEs also appear to clump slightly around the three major elliptical galaxies in Virgo: M87, M49 
and M86. Within a degree of each of these galaxies the velocity dispersion of dETs is essentially the same 
as the dispersion of all Virgo dETs. For the five dETs within 0.5° of M87, the dispersion is 428 km s -1 with 
a 14% change of being random, and an average velocity of 1097 km s^ 1 . The velocity of M87 is 1307 
km s" 1 , and thus there appears to be a strong clustering of a few dEs around it. The 18 dEs within 0.5° of 
M86 (with v = -244 km s _1 ) have an average velocity of 508±812 km s _1 , which has only a 0.1% change 
of being random. M49, located at the southern part of the Virgo cluster, is spatially a separate system, 
but has a similar velocity (997 km s _1 ) to that of the Virgo core region. We see some evidence for dET 
association with this elliptical galaxy and its surrounding cloud. Within 0.5° of M49 there are no dETs, but 
within 2° there are six dETs with an velocity of 1083 km s _1 and a low velocity dispersion of 291 km s _1 , 
with a 2% change of being random. This dispersion remains low (473 km s _1 ) for all dETs within 3° of 
M49. This has been known for some time for all galaxies in the M49 region (Binggeli et al. 1987). We 
conclude, in agreement with Ferguson (1992), that there is some evidence for bound dE companions to giant 
galaxies in Virgo, but it is only a statistically significant effect for M86. As is suggested from their spatial 
positions, most dEs are distributed in the general gravitational potential of Virgo and are not companions 
of larger galaxies. 

To determine if the velocity and projected locations of the various galaxy populations are statistically 
distinct, we perform Kolmogorov-Smirnov (K-S) tests on the positions and velocities of each individual core 
population against each of the other individual core populations. We do not perform these tests for the 
total galaxy populations, but only the core ones, to avoid biases from other Virgo galaxy clouds. These 
are two sided K-S tests that compare the distributions of velocities and positions of galaxies in different 
populations with each other (e.g., Press et al. 1992). Since K-S tests are based on statistics of two vectors, 
the positional test is limited to the projected radial distance from the center of the cluster. These tests are 
also somewhat constrained since if galaxies are in equilibrium their spatial positions and velocities should 
reflect the underlying mass distribution of the cluster. As such the interpretation of a lack of a correlation 
between two populations is not straight forward. On the other hand if two populations were accreted into 
Virgo at the same time, then we might expect their orbital and spatial positions to be correlated. 

The results of these tests are shown in Tables 5 and 6. Unfortunately, the results are largely inconclusive 
showing no statistically significant correlation between any two populations, but they are instructive in 
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several cases. We discuss individual K-S test for each Virgo galaxy population in the next section. 

3.2. Accretion and Orbit Signatures 

Figure 9 shows pie diagrams of slices through the total Virgo sample in both right ascension and 
declination. These figures have several features that illustrate the dynamical state of the Virgo cluster. The 
first is the clear 'finger of god' signature of the virializcd Virgo core. The population here is a mixture of 
all morphological types. 

Also in these figures are clear indications of galaxies still in the process of accreting into the Virgo core. 
These produce shell-like 'infall ring' regions, perpendicular to our line of sight in a pie diagram (Praton & 
Schneider 1994). In Figure 9 there are also several clumpy areas that contain galaxies, sometimes very close 
to the 'finger of god' feature. By far the most common galaxy to deviate from the 'finger' are the spirals 
which define these clumpy regions. The dETs (open circles) tend to follow both the spirals (X figures) and 
ellipticals (black squares); outer associations are dominated by spirals and dwarf ellipticals, but with a very 
high spiral/dET ratio. 

While velocities of spirals and ellipticals in Virgo are almost all known, only a very small fraction of 
dET velocities have been measured. Of these, most are for dETs in the densest regions, and mainly for the 
more easily observed dE,N types. It is therefore possible that dETs populate all regions of this diagram, 
especially in the outer associations dominated by later-type galaxies. Furthermore, velocity dispersions of 
various Virgo galaxies as a function of distance from the cluster core (Figure 10) indicates that dEs arc on 
complex orbits. This is also true for the spirals, SOs, dSOs and irregulars, indicating that these systems have 
a mix of orbits including highly elliptical or radial ones. 

3.3. Velocity Distributions 

If some galaxies in Virgo originate from an accreted component that underwent interactions with other 
cluster members, and if the velocities reflect only the gravitational potential well of Virgo at time of infall, 
then the various velocity dispersions of accreted populations should be somewhere between that of the old 
giant E component, assuming they trace the early potential (see §3.3.1), and currently infalling spiral and 
irregular galaxies. The evolution of the distribution of velocities of Virgo galaxy populations and their 
dispersions can be predicted by using reasonable assumptions about their total lifetimes in the cluster. This 
can be done by investigating the effects of two-body relaxation, equipartition of energy, and dynamical 
friction. 

To test these ideas we produced velocity histograms (Figures 4 & 5) and statistical data (Tables 5 & 
6) on the total (Figure 4) and core (Figure 5) Virgo galaxy populations, including: giant ellipticals (E), 
spirals (Sa-Sd), dwarf ellipticals (dE) nucleated dwarf ellipticals (dE,N), lenticulars (SO), dwarf lenticulars 
(dSO), and irregulars (dlrr, I, Sm, BCD). 

We fit Gaussian distributions to all of the core histograms (Figure 5) by using the Levenberg-Marquardt 
(LM) algorithm (Marquardt 1963). We use the results of these fits to test for kinematic characteristics 
of infalling versus bound galaxy populations. In particular, a Gaussian form is expected for a virializcd 
population. Thus finding how well a Gaussian can fit a velocity distribution is one method for determining 
if a population might be relaxed. The Levenberg-Marquardt method uses a non-linear least squares routine 
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to fit the number of galaxies (N) as a function of velocity (v) in the form: 



N = A + B xv + EIC xexp 



-2.7 x 



(v — center) \ " 
FWHM J 



(1) 



where A and B are constants, and always fit to zero, C is the amplitude of the Gaussian, while "center" 
and "FWHM" are the center and full width at half maximum for the fitted Gaussians. This summation 
is over the total number of Gaussians needed to fit the distributions. For cases where either one or two 
Gaussian components could be fit, we present in Figure 5 only single Gaussians. The centers of the peaks 
are also fitted by the LM algorithm, and not chosen a priori by us. 

The different velocity distributions for galaxy populations in the total Virgo sample can be seen 
in Figure 4. The ellipticals, spirals, lenticulars and irregulars have bimodal distributions. The second 
peak located near 2500 km s -1 is largely the result of projected clouds that are probably falling into the 
main cluster (Virgo A in Binggcli ct al. (1987) and Virgo I in dc Vaucouleurs & de Vaucouleurs (1973)). 
Interestingly, there is little evidence of this second 2500 km s _1 peak in the dET and dSO populations. This 
is an indication that nucleated dwarf ellipticals are abundant only in the well defined rich Virgo cluster 
proper, located near a velocity of 1000 km s -1 . The nucleated dwarf ellipticals and dSO population are 
associated with the dynamical "Virgo Core" or Virgo A cloud at V = 1064 km s _1 , although significant 
sub-structure is evident. 



3.3.1. Ellipticals 

The radial velocity distribution of ellipticals changes slightly between the total Virgo area and the Virgo 
core (Figures 4 and 5), with a 75% probably of association from a K-S test. The secondary distribution 
centered at 2500 km s -1 is almost completely absent in the Virgo core population, and the radial velocity 
distribution of the core ellipticals can be fit as a single Gaussian distribution with the best fit x 2 — 1-66, 
peak at 1093 km s -1 , and a velocity dispersion a = 462 km s _1 . 

Despite this, even in the core population, there are ellipticals with highly deviant velocities for their 
magnitudes (Figure 8). One galaxy at 2284 km s _1 , NGC 4168, is uncertain as a member of the Virgo 
cluster proper (VCC), and has an uncertain classification due to its probable low mass (Ho et al. 1997). 
Another high relative radial velocity galaxy is NGC 4473 with v = 2240 km s _1 . It is also doubtful that 
this galaxy is a proper member of the Virgo Core or even an elliptical. It contains disky isophotes, evidence 
of dust (Michard & Marchal 1994; van den Bosch et al. 1994) and its properties do not fit elliptical galaxy 
scaling relations (Davies et al. 1983). 

Excluding these two galaxies from the elliptical population we see a slight, but statistically insignificant, 
signature of equipartition of energy or effects of dynamical friction. The giant E galaxies appear by eye to 
have a decreasing dispersion with increasing mass, a possible signature of a process that has led to the giant 
ellipticals being more dynamically bound in the cluster. To test if this result is statistically significant, we 
carry out a Monte Carlo simulation of the velocities and magnitudes for these ellipticals. We characterize 
the mass- velocity segregation by measuring the difference between the velocity dispersion for the bright 
(My < —19) and faint (My > —19) ellipticals. The velocity dispersions of these two components are 299 
km s _1 and 509 km s , respectively, a 210 km s _1 difference. Using Monte Carlo techniques, we find a 2% 
chance that this difference is random, and that the observed difference is a 2.3ct event. It is therefore likely 
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that this magnitude-velocity difference is due to chance, and not a real physical effect. This signature is 
removed completely if we include the two rejected ellipticals. 

The ellipticals in the Virgo core region have the lowest total velocity dispersion a — 462 km s _1 of the 
major cluster populations, a Gaussian velocity distribution (Figure 5), and a centrally concentrated spatial 
distribution (Figure 6). These all suggest that ellipticals are the oldest and most relaxed component in 
Virgo. Distances derived from surface brightness fluctuations also suggest that most of the ellipticals in 
Virgo are gravitationally bound to the cluster (Neilsen & Tsvetanov 2000). Figure 10 shows the steadily 
decreasing value of the projected velocity dispersion as a function of projected distance from the dynamical 
center of Virgo for the elliptical population. Virgo ellipticals appear to be more dynamically relaxed in 
comparison to other populations. These effects are also seen for relaxed populations in N-body simulations 
of structure formation (e.g., Crone, Evrard & Richstone 1994; Thomas et al. 1998). 

3.3.2. Lenticular s (SOs) 

Depending on their nature, cluster lenticulars (SOs) may provide a critical clue to the evolution of 
galaxies in dense extragalactic environments (Dressier 1984). The major question to answer is whether or 
not this population has an origin more similar to the ellipticals or to the spirals. Being a 'transition-type', 
the lenticulars have properties of both ellipticals (old stellar populations, symmetric structure) and spirals 
(outer disks). 

The kinematics of the SO population in Virgo clearly differ from those of the ellipticals. The velocity 
dispersion of the SOs is 647 km s _1 , almost 200 km s _1 higher than that of the elliptical population, and 
similar to that of the dET population. Based on their radial velocities and positions, there is a 62% - 68% 
probability that the dEs originate from the same population as the SOs, although there is a lower 5% - 
39% probability for a velocity association with the dE,Ns. The SOs are clearly more extended spatially 
than are the ellipticals (Figure 6) and again resemble the spirals and dEs in this aspect. Their velocity 
distribution also contains a low velocity component near 300 km s _1 that is also seen in the dET, spiral 
and irregular populations. The overall velocity distribution for SOs is not fit very well by a Gaussian. The 
clumpy distribution of SO velocities indicates that these galaxies are not a relaxed Virgo component. 

Based on this we conclude that SOs are not associated with ellipticals, but formed later, potentially by 
a process similar to that by which the dEs form (cf. Mao & Mo 1998; Moore et al. 1999). This is consistent 
with observations of higher redshift z ~ 0.5 clusters, where the SO fraction is as much as three times lower 
than that seen in nearby clusters (Dressier et al. 1997), while the fraction of ellipticals in clusters is found 
to be the same at z < 1 (Dressier et al. 1997). High-resolution N-body simulations of cluster formation 
based on a CDM scenario also show that the properties of SOs cannot be reproduced from major mergers, 
while those of elliptical galaxies can (Okamoto & Nagashima 2001), suggesting that SOs formed later. 

3.3.3. Spirals and Irregulars 

The spiral galaxy population has a wide non-Gaussian velocity distribution with a total core velocity 
dispersion a ~ 776 km s _1 , over 300 km s _1 higher than that of the ellipticals (Table 4). The multi-peaked 
and broad velocity histogram (Figures 4 and 5) can be interpreted as a result of spiral-rich groups infalling 
into the cluster. Huchra (1985) characterized the velocity structure of the Virgo spirals as 'flat-top'. 
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However our lower errors allow finer binning, and there appear to be three components centered near 300 
km s _1 , 1100 km s~ x and 1900 km s" 1 . Similar peaks are seen in the dET, SO, and irregular populations. 
For an infall model, the high velocity components are spirals falling into Virgo on the side nearest the Local 
Group, while the low velocity components are spirals falling towards us on the opposite side. 

Further evidence that some Virgo spirals were recently accreted is shown by a pattern of lessening 
gas depletion within spirals at greater projected radii from the central cluster core, where the hot dense 
ICM could removed gas by ram-pressure (e.g., Haynes & Giovanelli 1986; Cayatte et al. 1990; Solanes et 
al. 2001). Spirals far from the core are mainly falling into the cluster, while those in the inner part of 
Virgo have begun to go through the core, and as a result of gas stripping by the ICM, are depleted of gas. 
Distances to Virgo galaxies derived from the Tully-Fishcr relation also reveals the expected infall velocity 
pattern: galaxies with blue shifts relative to the cluster are furthest away (Gavazzi et al. 1999). 

The spirals in Virgo are statistically similar in velocity space only to the irregulars (Table 5). Is it 
possible that irregulars were once satellites of present day Virgo spirals? During infall companions of spirals 
will be separated due to tidal forces. The potential of Virgo will induce, for pure linear motion and linear 
arrangement between, for example, a spiral of mass M sp and a companion, a maximum radial differential 
acceleration of 



a c iu — GM c i u 



((r + R)2) 



(2) 



The spiral and its companion will become dislodged when a c i u > GM sp /R 2 . In these equations, R is the 
distance between the infalling spiral and its companion, r is the distance from the leading edge of the group 
to the center of the cluster, M c i u and M sp are the masses of the cluster and infalling spiral, and G is the 
gravitational constant. This formula changes somewhat in a non-radial situation, since eq. (2) is only valid 
for a radial force. If we take the mass of Virgo as 10 M Q , concentrated within r, and an infalling spiral 
galaxy with mass 10 10 , then a galaxy separated by 0.5 Mpc from the host spiral will become dislodged 
within 2.5 Mpc of the cluster center. We can generalize this by using radial tides on an infalling galaxy of 
mass M ga i and a companion that will become dislodged when RGM c i u /r 3 = GM ga i/R 2 , which occurs at r 
= 10 Mpc. 

The irregulars^ in, or projected against, the core have characteristics very similar to the spirals and in 
some ways also to the dwarf ellipticals. The low and high velocity components are present in the irregular 
population (Figure 5) and the spatial distribution is very similar to that of the spirals (Figure 6). There is 
also a 92% probability that the irregulars originated from the same velocity population as the spirals, as 
shown by K-S tests (Table 5). 

From the triple peaked velocity distribution and the high velocity dispersion of 727 km s _1 , similar to 
that which we see for the spirals, it seems likely that the irregular population is infalling into the cluster 
along with the spirals (e.g., Gallagher & Hunter 1989). Their high gas content and star formation rates 
indicate these objects have not yet crossed the cluster core, where their extended gas will be rapidly 
stripped by the intracluster medium (cf. Lee et al. 2000; Mori & Burkert 2000). 



4 It is worth repeating for clarification that our irregular population includes Magellanic Spirals (Sm), as well as the classical 
dwarf irregulars and blue compact dwarfs. Any comparisons to other studies should keep this in mind. 
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3.3.4- Dwarf Ellipticals 

The overall dET population has a velocity distribution that is similar to that of the spirals but not to 
that of the ellipticals, with a velocity dispersion of 726 km s _1 , compared to the value of 776 km s _1 for 
the spirals, and 462 km s _1 for the ellipticals. This general trend has been suspected for some time and 
found by previous studies (e.g., Bothun & Mould 1988), but also has been disputed on the basis of smaller 
samples than the one used here (Binggeli, Tammann & Sandage 1987). The projected spatial distribution 
of the dETs, like that of the spirals, does not show as significant a decline in radius as does that of the 
ellipticals. To which of these two major populations are the dETs similar? The velocity distributions and 
dispersions of both dETs and pure dEs (shaded histograms in Figures 4 & 5) resemble more the spirals 
rather than the ellipticals. 

The distribution of dET velocities is not fit by a single Gaussian, another indication of dynamical 
evolution, youth and/or substructure. If Virgo dEs were in cquipartition with the giant ellipticals we 
would see evidence for mass-position segregation. In this case, the dEs would have a more extended spatial 
distribution, and a larger velocity dispersion than what is observed (§4.1.1). In fact, many dEs are in the 
richest part of Virgo, with fewer found in the lower density areas (Figure 6). Possibly some mass-position 
segregation is occurring, since the dEs are more spread out than the giant Es, but this is only a slight 
effect that could also be accounted for by dEs originating from accreted components not on radial orbits. 
The velocity distribution of the dEs is also clumpy, including evidence for a low- velocity peak (Figure 5). 
From these data, we conclude that the dETs taken as a single population are less dynamically evolved than 
the Virgo cluster core E populations, and are not relaxed. To further test this will require knowing the 
3-dimensional spatial distributions and velocity dispersion profiles of Virgo galaxies. 



4. Dynamical Models 

Modeling the dynamical evolution of galaxies in Virgo is critical for understanding these various 
kinematic characteristics and how they evolve. In this section, we investigate various ways the velocity 
distributions of galaxy populations can change over time. This includes determining what the velocity 
distribution would be if all galaxies have been in the cluster since its initial formation. Based on these 
calculations we conclude, consistent with our observational data, that the elliptical population is the oldest, 
and only relaxed sub-component of the cluster, while all the remaining populations have characteristics of 
later infall. 



4.1. Virialization and Relaxation 

The time for a galaxy in a cluster to relax can be approximated to an order of magnitude in several 
ways. One way is to use the two-body relaxation time (Chandrasekhar 1942; Spitzer & Hart 1971; Spitzer 
1987). For a self-gravitating population of identical mass bodies, the two-body relaxation time t r , defined 
as the time for a galaxy's orbit to become deflected by 1 radian, is given by 

1 I, . ^9 9 , 0.06N . . 

t r = a 3 /(47rG 2 m2n 9 lnA) « ^j^j X U (3) 
where a is the velocity of the galaxy, m g is the mass of the galaxy, n g is the number density of cluster 
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galaxies, A is the ratio of max to min impact parameters, N is the number of particles in the system, and 
id is the cluster crossing time. 

To compute the total relaxation time for the Virgo elliptical population we assume that early in Virgo's 
history only the ellipticals existed. Assuming a primordial average elliptical mass of 10 12 M , an initial 
number density <~ 20 galaxy/Mpc 3 , and a velocity dispersion of ~ 460 km s _1 (Table 2) gives a relaxation 
time of ~ 1.5 Gyr. The time for a system to completely relax is about 10 x t r , thus it would take ~ 15 Gyr 
for the elliptical population to complete this process. This is approximately a Hubble time, and thus it is 
unlikely that relaxation has completely occurred. 

The relaxation of the Es is however consistent with observations of galaxy populations in (z ~ 0.5) 
clusters where the 'core' ellipticals seem to be well in place (Dressier et al. 1997). Furthermore, the star 
formation history of moderate redshift clusters suggest that the stellar populations of ellipticals are mostly 
formed by z ~ 3 (Ellis et al. 1997). This formation must occur early to account for the mass-metallicity 
relationship found in nearby cluster ellipticals (Kauffmann 1996). 

This also agrees with results from numerical Cold Dark Matter models which show the first objects 
to relax in clusters are large elliptical galaxies (Diaferio et al. 2000). Simulations of cluster formation 
also indicate that the first galaxies formed in clusters should now be preferentially at the center (White & 
Springel 2000; Governato et al. 2001). The ellipticals are consistent with being these objects since they are 
the most centrally concentrated population in Virgo. 

For a dET population, similar to the present day density of several thousand dETs in a volume of a 
few Mpc 3 , the time for relaxation is > 10 13 yrs, assuming all the Virgo mass then was contained in dEs. 
Thus relaxation is not likely to have occurred for any dETs, old or young. 



4.1.1. Energy Equipartition 

We can crudely estimate features of dEs if they have achieved energy equipartition with the giant 
ellipticals. In this scenario the average energies of the elliptical and dwarf ellipticals should be the same. 
This would require that the average velocity of the dET population be (Mg/M^) 1 / 2 x Ve, where all 
quantities are averages. Taking the average mass of a Virgo elliptical and dwarf elliptical as 10 12 M Q , 
and 10 8 M Q respectively, and the observed velocity dispersion of 462 km s _1 for the ellipticals, the average 
velocity of the dEs should be ~ 46 x 10 3 km s _1 . This velocity is much larger than the observed 726 
km s _1 and is larger than the Virgo escape velocity. Dwarfs therefore would be escaping the cluster, and 
at the very least would have a very large velocity dispersion. We should also begin to see some mass 
segregation, with significant numbers of dETs in the outer parts of the cluster. Most dETs are however in 
the core region of Virgo, although with a large dispersion in distances (Young & Curric 1995; although see 
Binggeli & Jerjen 1998). 

The time scale for energy equipartition to occur is given by (e.g., Bertin 2000, eq. 7.10), 

ten = x tr, (4) 

m g 

where t r is the relaxation time (eq. 3) and mf is the mass of field particles. For the giant Es, t oq w t r = 1.5 
Gyr, although for dwarf ellipticals to reach equipartition with ellipticals t oq > 10 3 xt r ~ 10 7 Gyr. Based on 
this approximation the dETs and the other Virgo populations cannot be in equipartition of energy with the 
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cllipticals. This is verified by our observations (Figure 8) that shows no signs of equipartition in the dET 
or spiral population. However, the giant Es can be in equipartition with each other, if we assume that the 
Es were the first and only galaxies in the cluster several Gyrs ago, and all additional galaxies were accreted 
adiabatically. 

Any new dETs introduced into the cluster, as well as any old dETs, should show no signs of energy 
equipartition or full relaxation. Lacking these signatures therefore does not prove that Virgo dEs are a 
young cluster population. 



To understand how any new galaxy population introduced into Virgo will kinematically evolve, we 
consider the case of galaxies undergoing spherical accretion into a cluster at the turn- around radius, where 
the velocity of infall just cancels the expansion rate of the universe (Peebles 1980). 

In the special case where the density perturbation associated with a galaxy cluster has a constant 
density within the turnaround radius Rt, we can write the mass within any radius r, where r< Rt as 
M(<r) = (4/3)7r r 3 p. The infall velocity of a galaxy joining a cluster with mass M, and turnaround radius 
Rt scales as v 2 ~ px R„ ~ M/Rt- Thus, depending on the form of M/Rt, clusters with higher masses will 
contain accreted galaxies with higher radial velocities. This can generalized in the following way (e.g., Gott 
& Rees 1975). When the initial density enhancements in the universe can be represented by a power law 
after recombination, then 



where n > — 3 for hierarchical clustering. For a universe with fi = 1, the mass and size, Rt, of assembled 
areas are related by 



Unless M/Rt remains static, or decreases, which requires n<l, then galaxies accreted at later times when 
the cluster mass is higher, will have higher M/Rt ratios, and thus higher observed velocities than when the 
cluster contained less mass. This is consistent with spirals and irregulars being the most recently accreted 
component since they have the highest velocity dispersions. The rate of shell infall increases with higher 
51m, and for Qm = 1 all shells eventually fall into the cluster. The current Virgo mass within a radius of 
800 kpc is 1.8 - 3.3 xlO 14 M Q , based on fits to mass profiles (Girard 1999) while SBB99 find a somewhat 
lower value. In the past Virgo likely contained a lower total mass, and any galaxies accreted then will have 
lower observed velocities than those captured recently. 

Dynamical friction is also an important aspect to consider. We can compute the effective time scale for 
dynamical friction, in terms of the initial radius r*, velocity V and mass M, assuming a static potential, as 
(Binney & Tremaine 1987): 



4.2. The Evolution of Velocity Structures and Evidence for Infall 




(5) 




(6) 
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For the ellipticals, tf r ~ 3 Gyr. This is short enough for dynamical friction to have an effect on the kinematic 
structure. This is especially true if ellipticals have existed in clusters for many Gyrs, as is suggested by 
observations (e.g., Ellis et al. 1997) and theory (Kauffmann 1996). For all other lower mass populations, 
including the dEs, dynamical friction is not an efficient process with t fr <~ 10 3 - 10 4 Gyrs. 

We conclude from this argument, and the previous calculations in §4.1, that velocities of low mass 
(< 10 10 M Q ) accreted galaxies, and hence their velocity dispersions, should not change very much over a 
Hubble time due to dynamical friction, relaxation or energy equipartition. If dEs are a new population 
that resulted from infall in the last several Gyr, then they should retain many of their original orbital 
characteristics. Two possibilities then exist for explaining a separation between giant ellipticals and the 
dEs. The ellipticals and dEs could have formed together and later the distribution of ellipticals modified to 
become more centrally concentrated due to dynamical evolution. Alternatively, the properties of the dEs 
may reflect later infall into the cluster. 

We can use the fact that the ellipticals are relaxed while the other populations should retain their 
original orbital features to argue that most spirals and irregulars have recently been accreted, while the 
dET and SO populations primarily originated from galaxies that fell in after the cluster core formed. For 
a virialized population the kinetic (T) and gravitational potential energy (U) should be related by I T I = 
1/2 |UI, while for an accreted population U + T = and thus |T| = |U|. Assuming that on average 
the potential of a population that fell into Virgo is equal to the potential of a virialized population, we can 
relate the velocity dispersions of a virialized and infalling population, assuming isotropic orbits, as (see also 
Colless & Dunn 1996) a in f a n = \/2 x a v i r = 1.41 X(j vir . For the Virgo populations the ratios of <T pop /<JE, 
where pop is any one of dEs, SOs, irregulars, or spirals, are ct^et /&e = 1-57 or 1.51 if we include dSOs, 
cso/ce = 1-4, crirr/fE = 1-57, crs p i ra i/(TE = 1-65. These values are all very close to that expected for the 
ratio between infalling and virialized populations, and are not obviously accounted for by the dynamical 
evolution models. We therefore favor the delayed infall model for the source of cluster dEs. If infall is 
occurring, orbital characteristics should also correlate with galaxy types; for example, star forming galaxies 
should have unique orbital properties (Mahdavi et al. 1999; Balogh, Navarro, & Morris 2000), including 
large velocity dispersions which we see (Figure 10). 



4.3. High-Speed Interactions 

Although neither dynamical friction nor relaxation induce major effects on the evolution of kinematics 
of low-mass cluster galaxy, high-speed or impulse interactions are potentially an important method for 
altering the internal properties of infalling galaxies. The interactions between galaxies in Virgo, with a total 
velocity dispersion of <~ 700 km s _1 , are nearly all at high-speed, at relative velocities much higher than 
individual internal galaxy velocity dispersions. Thus, it is instructive to investigate how these interactions 
can affect the properties of cluster galaxies. As we show here, high-speed interactions, while important, are 
rather inefficient at destroying low-mass galaxies, but are likely to significantly modify their structures. 

We can approximate the strengths of high-speed interactions by using the impulse approximation (cf. 
Spitzer 1958; Binney & Tremaine 1987). The energy imparted into a galaxy during a non- interpenetrating 
impulsive interaction can be computed in the tidal approximation by (Spitzer 1958), 
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where Mi and M 2 are the masses of two systems undergoing the impulse, b is the impact parameter, V is 
the initial relative velocity of the two systems and r 2 is the mean square radius of the perturbed system. 
This energy typically goes into the individual stars that can then be liberated from their host galaxies 
(Gallagher & Ostriker 1972; Aguilar & White 1985). 

An important question to ask is if this internal energy addition will be enough to destroy individual 
galaxies. Figure 11 shows the ratio of AE to the internal energy, defined as 1/2 Mi a 2 , of a galaxy 
interacting with a M = 10 12 M system, as a function of the impact parameter, b from eq. (8). Both lines 
represent systems with internal velocity dispersions a = 50 km s _1 and relative velocity V = 1000 km s _1 , 
with the solid line a system with radius r = 5 kpc (a spiral) and the dashed line with radius r = 2 kpc (a 
dwarf). The ratio of the increase in internal energy from an impulsive interaction to the internal binding 
energy of a galaxy with velocity dispersion a is given by, 



AE 8 G 2 M 2 2 



(9) 
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Figure 11 shows the relative unimportance of this energy increase for destroying galaxies undergoing a 
single impulse. Typically only dwarf systems that come within 10 kpc of large ellipticals with mass ~ 10 12 
M will receive enough of an increase in internal energy to disrupt the system in one encounter. These 
distances are rarely, if ever, reached in a cluster of galaxies such as Virgo. In addition, at this close distance 
the impulse approximation breaks down, and these equations no longer hold. The impulse energy from 
high-speed interactions is thus unlikely to destroy the large number of dwarf ellipticals in Virgo at least 
on short time scales. This will, however, potentially have an important effect on each galaxy's internal 
evolution (see §5.3). 



5. Discussion 

In the previous section we demonstrated how dETs, along with the Virgo spirals, irregulars and SOs, all 
have kinematic signatures of past accretion. In this section we discuss several possible origins of the Virgo 
dETs. This includes investigating whether or not dETs could be accreted field dwarfs, or if they formed 
from some dynamical process induced by the cluster. 



5.1. Infalling Field and Group Dwarfs 

Could a significant fraction of Virgo dETs originate from field or group dEs accreted into Virgo? Is it 
possible that some or all the dETs in Virgo were once part of something like the Local Group (LG) that 
fell into Virgo, with the dETs dislodging from the spirals? The answer depends on the number of dwarfs 
in groups and in the field, an uncertain quantity (e.g., Ellis et al. 1996; Trentham 1998). We do however 
have some limits on the number of field dwarfs, and particularly the number of dwarf galaxies around 
spirals, such as those in the Local Group. The number of dwarf galaxies per giant is very high in clusters, 
much higher than in the field (e.g., Binggcli ct al. 1990; Seeker & Harris 1996). If all the spirals in Virgo 
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originatcd in systems similar to the LG then in the magnitude limit of the VCC (By < 18; Mb < —13.3 
assuming m-M = 31.3) there are only five LG dwarf spheroidals that could be be detected. These are NGC 
147, NGC 185, NGC 205, the Fornax dSph, and the Sagittarius dSph. We consider M32 a stripped giant 
elliptical, and not a dwarf elliptical, for these purposes, although including it as a dE does not change the 
following results. 

There are only three bright LG spiral galaxies and they are all detectable with the same magnitude 
limits: The Milky Way, M31 and M33. This gives a dwarf elliptical to giant spiral ratio of 5/3 = 1.67. In 
the VCC there are 1277 members, with 574 more possible members, most of which are dwarf ellipticals. 
There are, however, only slightly over 100 spiral galaxies in Virgo (<~ 160 if we include SOs), which would 
correspond to less than 300 detectable dwarf ellipticals if the dET/gS ratio were the same as in the LG. 
Similarly, if the 40 Virgo cluster ellipticals originated from mergers of spirals, they would also come from an 
initial population of ~T00 spirals, and we would again expect to find only a few hundred dE companions. 
The VCC has well over 1000 dEs at these limits, and thus a simple scenario where the Virgo dEs are old 
members of galaxy groups dislodged from their spiral parents under predicts the number of dEs by over a 
factor of 3. There is also no sign of excess dEs associated with currently infalling spirals. Since the ratio of 
E/SOs to dEs is lower than that which see in groups we conclude that dETs in clusters cannot be accounted 
for by accretion of field dEs, or dwarfs attached to infalling spirals that later become removed. 

5.2. Gas Stripping: dl -> dE 

One of the historically more popular methods of creating a dwarf elliptical is by stripping gas from 
a star-forming dwarf irregular (Lin & Faber 1981). This type of evolution is possible in Virgo since the 
dls, along with the spirals, have an infalling signature and thus it is possible that a large number of dls 
were accreted in the past. We do know that the large number of dETs cannot be accounted for by the 
relatively small number of present day dls (Gallagher & Hunter 1989), thus a much higher past dl accretion 
rate would be required. We therefore conclude that the vast majority of the dETs in this study, that are 
by necessity bright, did not evolve from irregulars. However, because of the huge population of cluster 
dETs and their large range in sizes and magnitudes, it seems possible that some of them, particularly the 
low mass ones, evolved from irregulars. Irregulars have low masses and it is possible that after rapid gas 
stripping processes such as impulse encounters would produce a very small and faint dwarf elliptical after 
many Gyr. It would be very useful to extend searches for an extremely faint population of dEs that should 
exist, as products of the evolution of dls that have been accreted into clusters over many epochs along with 
their spiral companions (cf. Impey, Bothun & Malin 1988). Some of these have potentially been found in 
Virgo (Caldwell & Armandroff 2000). 

5.3. Galaxy Harassment: S — > dE 

Galaxy harassment (Moore et al. 1996; Moore et al. 1998) is the process whereby cluster galaxies are 
stripped of their interstellar medium and become dynamically 'heated' by high speed interactions with other 
cluster galaxies and the cluster's gravitational potential. During these encounters the internal potential 
energy of an infalling galaxy increases creating a situation where the galaxy is no longer in equilibrium. To 
re-establish equilibrium the galaxy increases in size and loses its most energetic stars. Because of this, the 
galaxy's density decreases over time and after several orbits a spiral can morphologically transform into a 
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dwarf elliptical (Moore et al. 1998). The lost mass necessarily becomes part of the intracluster material. 

Galaxy harassment can explain a number of observations. Moderate redshift clusters at z ~ 0.8 have 
a high fraction of spirals and other star-forming galaxies (e.g., Butcher & Ocmlcr 1978; Butcher & Ocmlcr 
1984; Dressier et al. 1994; Oemler, Dressier & Butcher 1997) that are not seen in the same frequency in 
nearby clusters. These blue galaxies typically do not have identifiable merging or interacting companions, 
and therefore could be products of high-speed interactions, or galaxy interactions with the cluster potential 
inducing star formation. Features in nearby clusters, such as debris arcs (Mobasher & Trentham 1998; 
Gregg & West 1998), and distorted galaxies without obvious companions (Conselice & Gallagher 1999), are 
consistent with the harassment ideas of cluster galaxy evolution. The intracluster light in nearby clusters 
has similar colors to dwarf ellipticals, suggesting that material similar to dwarfs, perhaps the remnant 
material from stripping, is the source of this light (e.g., Seeker et al. 1997). Stripping is further confirmed 
by finding red giant branch stars (Ferguson, Tanvir, & von Hippel 1998) and planetary nebula (e.g., Mendez 
et al. 1997) in the Virgo cluster intragalactic regions. 

The kinematic results presented here also support the galaxy harassment scenario for a dET origin. As 
shown in §4, the velocity characteristics of low-mass (< 10 10 M Q ) infalling galaxies should change little 
over several Gyrs and thus their velocity and spatial distributions should reflect the mass profile of Virgo 
at the time of accretion (Carlberg et al. 1997). Galaxies accreted by more massive clusters will have higher 
observed velocities than those galaxies accreted by a lower mass cluster of the same size. If galaxies in 
Virgo were distributed in a similar way, those accreted early will have existed in the cluster longer, and 
thus should be more dynamically stripped by high-speed interactions than those accreted recently. In this 
scenario accreted galaxies with the lowest masses should have the lowest velocity dispersions. After the 
giant ellipticals, the dSOs and dEs have the lowest velocity dispersions (Table 4) with a combined a = 621 
km s _1 . The SOs and dE,Ns, generally more massive than the dSOs and dEs, have a combined a ~ 705 
km s _1 , over 80 km s _1 higher than that of the lower mass dEs and dSOs. This difference and the difference 
in spatial distributions is potentially reflecting the different masses of Virgo when each galaxy type entered 
the cluster. This is also consistent with the spirals and irregulars infalling now, since these galaxies have a 
combined a ~ 756 km s _1 , higher than any of the other galaxy types. 

This, however, is circumstantial, since the cluster mass given by each population depends upon their 
spatial and velocity distributions. This can be investigated in detail by using the Jeans equation that 
relates orbital and spatial distributions of particles in a spherical system to that system's mass profile M(r) 
(Binney and Tremaine 1987; Carlberg et al. 1997), 



where v(v) is the radial density profile parameter, v 2 is the radial velocity dispersion and (3 is the anisotropy 
parameter, [5 = 1 — Vg/v 2 . It is impossible to derive the value of /? without knowing the mass profile of 
Virgo, or obtain the mass profile with out an assumption about (3. The luminosity density v, and total 
velocity dispersion, v 2 are real space quantities, while we observe the projected luminosity density N and 
projected velocity dispersion a 2 . Figure 12 shows the projected spatial distribution of various populations, 
and Table 7 list the projected slopes for both the velocity and spatial distributions. It is possible to convert 
these slopes into three dimensional quantities by using the Abel transformation (e.g., Binney & Mamon 
1982), although we do not have enough information to do this. Additionally, only the dETs have a projected 
spatial distribution slope that can be measured with certainly out to a few degrees in radius. 




(10) 
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Other recent evidence of S — > dE evolution includes luminosity profiles of dEs and dE,Ns (Stiavelli 
et al. 2001), demonstrating that dETs have inner profiles more similar to spiral bulges than those of 
giant ellipticals. Many Virgo dETs also contain disk-like, as opposed to elliptical-like profiles (Ryden et 
al. 1999) and there are examples of Virgo dETs with faint spiral structure, including IC 3328 (Jerjen, 
Kalnajs & Binggeli 2000). There are also Virgo dETs with dust and gas features commonly associated with 
late- type star forming galaxies (Elmegreen et al. 2000). These observations all suggest that dETs could be 
descendents of star forming disk galaxies. 

How could so many cluster dEs (VCC) originate from spirals? The rate of infalling galaxies into clusters 
is predicted to peak at z <~ 0.8 in a CDM cosmology with Qm = 1, Hq = 50 km s _1 Mpc -1 (Kauffmann 

1995) . This is the same redshift where the Butchcr-Ocmlcr effect is found and where the infall rates into 
clusters is observed to be high (Ellingson et al. 2001). If dEs originate from spirals, then a significant 
number of them should have been introduced into the cluster at approximately this redshift. Since dEs are 
galaxies with infall signatures, including a non-Gaussian velocity distribution, that cannot be intrinsic to 
an old cluster population, or to accreted field dEs, they likely originate from this epoch. They are therefore 
possibly the modern descendents of harassed Butcher-Oemler galaxies. 

5.3.1. dEs vs. dE,Ns 

The velocity characteristics of dEs and dE,Ns are not significantly different (pure dEs are shaded in 
Figures 4 & 5). This suggests that these two populations have similar origins and that dE,Ns are not 
scaled down giant ellipticals. The one possible difference between dE and dE,Ns is their spatial distribution 
(Figure 6), with the dE,Ns appearing more centrally concentrated in Virgo (Ferguson & Sandage 1989). 
Analysis of VCC positions for the dEs and dE,Ns shows that although the average distance from M87 is 
different by 0.31° for the dEs (2.25±1.73°) and the dE,Ns (1.94±1.64° ), this difference is not significant 
when we consider the large ~ 1.7° dispersions, and when errors are taken into account. 

Ferguson & Sandage (1989) found a luminosity dependence in clustering for Virgo dEs, such that 
dE,Ns and the faint non-nucleated dEs are centrally concentrated, but the bright non-nucleated dEs are 
not. There are several possibilities for explaining these trends. The nuclei of dE,Ns could be induced from 
gravitational effects due to the cluster core. This can occur in two ways. One could be a gravitational effect 
similar to the one producing the large fraction of barred spirals found in the Virgo core region (Andersen 

1996) , where tides could induce an increase in stellar density at the cores of dwarf galaxies. On the other 
hand the gravitational potential could also stabilize dEs in the central regions of Virgo, allowing globular 
clusters to accrete to the center by dynamical friction, forming nuclei (Hernandez & Gilmore 1998; Oh & 
Lin 2000; Lotz et al. 2001). In addition, about 20% of Virgo dEs have nuclei offset from their centers 
(Binggeli, Barazza & Jerjen 2000). This is a possible indication that external gravitational fields are driving 
the nuclei to oscillate about galaxy centers, and thus playing a role in dE nuclei evolution. 

Another possibility is that nucleation is produced by star formation induced in the centers of these 
galaxies after they are accreted into the cluster. Any spirals that orbit close to the center of the cluster 
will experience gravitational torques which could induce gas to inflow into a dE,N progenitor's center 
producing bursts of star formation. Nucleated dEs also have a higher specific globular cluster frequency 
than non- nucleated dEs (Miller et al. 1998), an observation explainable by increased star formation in dE,N 
progenitors that orbit close to the Virgo core. 

These observations suggest that the harassment process of transforming spirals into dEs is likely 
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occurring, but probably is not the only source of these galaxies. Since most galaxies are in groups, and 
these fall into clusters, it seems certain that some dEs in Virgo were formed outside the cluster in galaxy 
groups. Our observations similarly do not rule out the possibility that some dEs in the core of Virgo were 
present when the cluster Es formed. 

6. Summary and Discussion 

In this paper we derive the following conclusions based on the kinematic properties of Virgo cluster 
galaxies: 

I. Elliptical galaxies in the Virgo core region form a relaxed, or nearly relaxed system. This is demonstrated 
by the Gaussian velocity distribution of the ellipticals, and their centrally concentrated spatial distribution. 
No other galaxy population has characteristics of relaxation. This includes the SOs, dEs, dE,Ns, dSOs, as 
well as the spiral and irregular galaxies. Relaxation and dynamical friction are shown to be inefficient 
for significantly decreasing the velocity dispersion of galaxies with masses lower than 10 10 M Q , although 
these are potentially important for the ellipticals. Mailing galaxies should therefore retain their kinematic 
signatures even after many cluster crossings. 

II. All populations, except for the ellipticals, have similar velocity dispersions (Table 4), with an average 
Cpop/<7E = 1.49, close to the expected ratio of velocity dispersions for a marginally bound and virializcd 
population. The dET, SO, spiral and irregular galaxies have spatial distributions that are progressively less 
centrally concentrated and have increasingly non-Gaussian velocity distributions. We therefore conclude 
that most galaxy populations in Virgo, besides the giant ellipticals, have been accreted into the cluster after 
the giant Es were in place. 

III. We investigate several possible origins for the dEs in Virgo including: dislodging of dEs companions 
connected to infalling spiral hosts, transformation of infalling dls into dEs, and stripping of material from 
infalling galaxies through processes such as harassment. We favor the last scenario with its implication that 
accreted spiral galaxies in the past were transformed into modern cluster dETs. This is suggested in part: 
the least massive galaxies, and those predicted to be the most stripped (dSOs and dETs) have the lowest 
velocity dispersion (621 km s -1 ) of any population besides the large ellipticals. The spirals and irregulars, 
the newest additions to the cluster, have a combined a = 756 km s _1 , the highest for any population, 
potentially reflecting the higher mass of Virgo when the dE progenitors were accreted. To confirm this will 
require a better understanding of the velocity anisotropy of the various galaxy populations and their three 
dimensional spatial distributions. 

These results have several implications. First, most dETs in clusters are not primordial cluster 
galaxies, or galaxies that have existed since the large ellipticals were in place. Since dETs are common in 
clusters, have dynamical properties that suggest they were accreted by Virgo several Gyrs ago, and were 
more massive in the past, they are potentially the remnants of Butcher-Oemler galaxies. This raises the 
question of what happened to the original Cold Dark Matter objects? Did they all merge early to form 
the giant elliptical population, or is there a large population of objects with masses and sizes similar to 
dEs, but with low surface brightness, that remain undetected (Cen 2001), or was the formation of low-mass 
halos suppressed (Bullock, Kravtsov & Weinberg 2001)? Either solution for explaining away the lack 
of detected old low-mass cluster objects has deep implications for observational cosmology and deserves 
further attention. 
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This study suggests that many cluster dEs are fundamentally different from dwarf ellipticals, or dwarf 
spheroidals, found in groups, although they could both be produced by harassment (Mayor et al. 2001). 
Establishing if this is true will require more observational evidence to show conclusively, but there may be 
two ways to form dEs - a simple galaxy collapse (e.g., Dekel & Silk 1986), and stripped down galaxies as 
suggested here. This paper also demonstrates how galaxy evolution can affect cosmological measurements. 
Galaxy clusters would be brighter if harassment did not occur, since galaxies are stripped of their stellar 
material. As a result this can change our estimate of cosmological biasing in dense regions. 

Future papers in this series will explore the properties of the stellar populations in cluster dEs 
themselves to learn more about their histories, and to identify exactly what their ancestors were. If cluster 
dEs arc "primordial" cluster members then their stars should all have roughly the same age, with no 
star formation after their initial creation, due to ram-pressure effects. If cluster dEs form by some other 
process, such as harassment of infalling spirals that were once BO galaxies, then there should be significant 
intermediate or young stellar populations. The stellar population ages of dEs should also correlate with 
their orbital kinematics. 
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TABLE 1 

Velocity Standards 



Name 


R.A. (J2000) 


Dec. (J2000) 


Spec. Type 


Velocity km s 1 


HD 29139 


04:35:55.2 


+16:30:33.5 


K5III 


54.1 


HD 62509 


07:45:19.0 


+28:01:34.3 


KOIIIb 


3.3 


HD 84441 


09:45:51.1 


+23:46:27.3 


GUI 


4.8 


HD 135722 


15:15:30.2 


+33:18:53.4 


G8III 


-12.3 


HD 136512 


15:20:08.6 


+29:36:58.4 


KOIII 


-53.1 


HD 137422 


15:20:43.7 


+71:50:02.5 


A3III 


-3.9 
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TABLE 2 

New Virgo Cluster Radial Velocities 





"R 4 (+2000 C\) 

xL.iT.. \o ^yjyjyj t\j j 


Dec H2000 01 


j_ype 


lVf a cmitn Hp 

1VJ.CX tl±l± (j UU.C 


v f 1c TY1 / 1 

V I / O I 


T? -vfl lnp^ 


VCC 292 


12:18:31.4 


+05:51:09 


dEl,N 


17.0 


2232+66 


3.2 


VCC 300 


12:18:42.2 


+05:40:03 


dE5 


18.1 


98541+61 


3.0 


VCC 303 


12:18:42.8 


+05:47:27 


dE2?,N 


15.8 


2846+43 


5.6 


VCC 314 


12:18:57.3 


+06:11:05 


SBc? 


14.9 


5724+40 


3.9 


VCC 315 


12:19:00.5 


+06:05:40 


SBa(s) 


15.0 


1525+31 


7.9 


VCC 321 


12:19:05.6 


+05:44:27 


dE6 


18.0 


1227+73 


3.3 


VCC 323 


12:19:06.4 


+05:43:33 


Sab 


14.9 


3002+39 


5.2 


VCC 336 


12:19:17.6 


+05:52:39 


dE,N 


16.2 


2990+42 


4.7 


VCC 338 


12:19:20.2 


+05:28:00 


dE3,N 


16.5 


2802+53 


4.3 


VCC 362 


12:19:42.2 


+05:32:17 


SB(s)a 


14.6 


1617+40 


5.6 


VCC 366 


12:19:45.3 


+06:00:21 


SO 


15.0 


2158+41 


7.3 


VCC 390 


12:20:04.4 


+05:24:57 


dE3 


16.9 


2479+38 


5.4 


VCC 684 


12:23:57.7 


+12:53:13 


dE,N 


16.2 


537+58 


4.7 


VCC 719 


12:24:19.0 


+12:54:47 


dE2 


18.5 


32738+59 


3.1 


VCC 723 


12:24:22.1 


+13:01:37 


dSO? 


15.0 


-111+25 


9.6 


VCC 753 


12:24:51.6 


+13:06:40 


dE,N 


16.3 


31166+57 


3.8 


VCC 762 


12:25:02.8 


+07:30:23 


dE3,N 


15.3 


1200+25 


6.7 


VCC 765 


12:25:03.5 


+13:14:40 


dEl,N 


16.0 


854+36 


6.8 


VCC 781 


12:25:15.0 


+12:42:52 


dS0,N 


14.8 


-188+23 


13.2 


VCC 817 


12:25:37.7 


+15:50:06 


dEl 


15.0 


1086+57 


3.4 


VCC 838 


12:25:46.9 


+12:45:37 


dE2 


17.5 


88564+72 


3.8 


VCC 854 


12:25:55.6 


+12:46:11 


dE8,N 


17.3 


684+75 


3.8 


VCC 871 


12:26:05.6 


+12:33:34 


dE,N 


15.4 


1636+76 


3.1 


VCC 872 


12:26:00.7 


+12:51:40 


dEO.N 


17.0 


1183+50 


3.7 


VCC 896 


12:26:22.5 


+12:47:00 


dE3,N 


17.8 


2114+63 


2.7 


VCC 916 


12:26:33.3 


+12:44:24 


dE,N 


16.2 


1115+32 


6.9 


VCC 928 


12:26:39.6 


+12:30:48 


dE? 


16.2 


-254+54 


4.9 


VCC 967 


12:27:04.0 


+12:51:54 


dE4 


18.7 


1135+51 


3.1 


VCC 977 


12:27:11.1 


+12:02:18 


dE4,N 


17.9 


102+45 


4.2 


VCC 996 


12:27:21.2 


+13:06:36 


dE5 


18.4 


-28+32 


7.7 


VCC 1027 


12:27:38.0 


+12:52:48 


dE0,N 


18.1 


77+46 


5.2 


VCC 1040 


12:27:44.5 


+12:58:55 


dE3,N 


17.5 


-143+63 


5.2 


VCC 1042 


12:27:45.8 


+12:52:19 


dE4 


18.6 


-62+13 


19.4 


VCC 1237 


12:29:51.2 


+13:52:03 




15.5 


-255+72 


2.9 


VCC 1351 


12:31:17.5 


+13:49:42 


dE4 


16.0 


187+48 


6.2 


VCC 1870 


12:41:15.3 


+11:17:55 


dE6 


15.8 


1617+27 


8.7 


VCC 1881 


12:41:29.5 


+10:45:46 


dE3;,N 


16.2 


138+56 


3.9 


VCC 1891 


12:41:48.9 


+11:11:29 


dE4,N 


16.7 


1016+43 


4.6 


VCC 1919 


12:42:19.0 


+10:34:04 


dE0,N 


17.0 


1869+54 


3.2 


VCC 1927 


12:42:35.9 


+10:33:55 


Sc 


15.0 


20199+14 


9.6 


VCC 1948 


12:42:58.0 


+10:40:55 


dE3 


15.1 


1581+22 


9.8 


VCC 1958 


12:43:10.0 


+11:02:12 


dE2,N 


17.0 


1049+30 


3.5 


VCC 1971 


12:43:31.0 


+11:02:50 


dE3: 


16.6 


1376+34 


6.1 



a Morphological types are from the estimates given by Binggeli, Sandage & Tammann (1985) 
b See Section 2 for explination 
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TABLE 3 

Virgo Dwarf Elliptical Galaxies with Radial Velocities 



Name 


ID 


R.A.(J2000.0J 


Dec. (J2000.0J 


Typo A 


Magnitude 


v(km/s) 


Source 


VCC 9 


TO 301 Q 


1 9-0Q-99 


tio. oy . 00 


dEl N 


13 93 


1 80,1 4- zlQ 

!OU4l4y 




VCC 31 




1Z. 1U. OO. I 


4-00- 1 ^-07 
Tuy. xo. u 1 


dE 


14 87 


991 c j_/in 

Z Z X OX4U 


2 


VL/Vj IDO 




X Z . 1 . i>4 . O 


4-1 zl-01 0f\ 


dE2 


17 


R89-I-1 
ooz zn n.) 




\ \^\^ 1 i o 


to '}nsii 


1Z. lO . U'.l . -k 


9. A 1 - OR 
-|-lZ.41.ZO 






1 IDoztOO 


1 


vP p 900 
V ^ ^ zuu 




i o. i ft . o o n 
1Z. ID. OO. y 


4-1 '3 ■ 01 ■ ^ft 
-\- ±oA)±. OD 


dE2'N 


Izt'fiQ 






VP P 91ft 
VIjIj Z10 


TO '}OQ'7 


19.17.n1 

1Z. 1 1 . Ul . Z 


_X.nO- 9/1 ■ '39 
-|-Uy . Z4t. OZ 


Q Hi O , p , 1 N 


a a 


1 '39 n, 4_ 70 
lOZoit ( y 




VP P 007 
Vl(^ ZZ 1 


T TO O 

1. ^0±4T 


12" IT" 14 4 


iria. r: ft. 09 
-|-Uo. OD. OZ 




on 


1 9004-01 

izyuzt y 1 




vp p oin 
V \s O ZOU 


IC 3101 


12" IT" 19 5 


4-11- ^.ft ■ '39 
-4- 1 1 . OD. OZ 




r o 


1 '3 00 4_7'J 

ioyo± ( 




VP 9Q9 

v zyz 






inr; . c 1 . no 

~l-uo. 01 . uy 


dEl'N 


17 


Z Z O Z It D D 




V ouo 




19-1 8-zL9 8 

X Z . IO.IZ.O 


4-0^-47-97 
A^ U . 1 1 . z ( 


dE2^ N 


15 8 


98zLfi4-zL ,: t 

Z040I40 




v oiy 


IC 3142 


12" 19"01 7 


TlO. OO. OD 


rlPO M 
t.L Xj'^ , 


16 2 


-206 




VPP Q91 

V V\J Oil 




j_z. ly . uo. 


4-0 c .-zlzL-97 

TU0.44. Z ( 


dE6 


18 


1 9974-70 
1ZZ i Az 1 O 




VPP 




12" 19" 17 6 


1 nr.. c.9- ^0 

tuo. oz. oy 




16 2 


9Q004-zL9 
zyyux4z 




vpp 'vm 

V^Vj ooo 




19-1 0-90 9 
1 z . 1 y . zu . z 


4-0^-98-00 
tuo. zo. uu 


dE3 N 


16 5 


zouz rjzijo 




vpp ■s.qo 




1 9-90-OzL A 
1 Z . Z U . u^± . ^± 


4_0 c .-9zL- ^.7 

TUO. Z4t. O 1 


dE3 


16 9 


9/1704- 'iR 
zi * ynzoo 


3 


VPP '3.07 




1 9-90-1 9 9 
J. Z . ZU . X z . z 


U D . O 1 . Z O 


dE5 N 


15 


9/1 nr. _|_ ft n 
z4yon:uu 




vpp zio^ 

V «w. 4UO 




1 9-90-1 8 1 

X Z . Z U . X O . X 


4-0^- ^.7 
tuo. oy . i 


dEO 


20 


2097 




VPP A ^7 


TTOO 7'-iQQA 


1 9-90-zLS 1 

IZ.ZU.IO. X 


+ 17"29" 16 


dE5 N 


14 54 


1 A7A-\-Afi 


2 


VPP r izL9 

V 04Z 


TO ,; t91 8 


1 9-99-1 Q R 
x z . z z . ly.o 


tuu. 00. oy 


dEO N 


15 00 


1 n r .74- ,; i0 
iuo 1 znoy 


r ( 


VPP ^.A^ 


TTOO 7A'Zi\ 


1 9-99-1 Q 

X Z . Z Z . ly.O 


4-1 A-A r \- 'iQ 
xi.io. oy 


dE5 


14 77 


861i58 


2 


VPP r./lc 

V \-S ^_/ 04:0 


TO ^H^tA 


1 9-99- 1Q K 
xz. zz. xy . D 


4-1 ^- /I /I ■ 09 
-|-10.44. UZ 






1 ioy it'iy 




VPP ^./l ft 




1 9-99-91 ft 
XZ . ZZ . Z X . 


11 n. '3ft. r .7 
-|-1U.0D.U ( 




15.7 


90ft74- 1 O/l 
ZUD ( it 1U4 




vpp fins 


MOO A^OO 


1Z. ZO. Ul . 1 


-|~ 15: 54: 21 


dE4, JN 


il o!i 


1 on '3 4- 1 OO 




VPP o /I 
V \^ D041 


MOO /I '}9Si 


1 9 . 9 Q . 1 Q Q 

xz. zo. xy . y 


-|~ 15:49: 13 


oIpo'w 


■ 


lyy-jzZ'i: 




VPP 

V \^ O 084 




1 9- 9'? ■ 7 
XZ. ZO. O i . i 


-(-12:53: 13 






537-jz58 




VPP 7/1 r. 


MOO AQfift 
In^jO 4TODD 


1 9-9/1 ■ A 7 O 
1Z.Z4.4 ( . U 


4_07- 91 ■ 1 o 

-|-u 1 . zx . xu 


dEf'^N 


14 67 


1 O Q ft 4- '3 

izoDzt oy 




VCC 750 




12:24:49.6 


+06:45:33 


dE5,N 


15.0 


1083+38 


5 


VCC 762 




12:25:02.8 


+07:30:23 


dE3,N 


15.3 


1200+25 


3 


VCC 765 




12:25:03.5 


+ 13:14:40 


dEl,N 


16.0 


854+36 


3 


VCC 786 


IC 3305 


12:25:14.5 


+ 11:50:57 


dE7,N 


15.11 


2388+30 


2 


VCC 790 




12:25:17.5 


+ 14:10:21 


dEl,N 


16.4 


573±54 


8 


VCC 797 




12:25:24.2 


+ 18:08:29 


dE3,N 


17.0 


773 ±8 


9 


VCC 810 




12:25:33.6 


+ 13:13:37 


dE0,N 


16.5 


-340+50 


1 


VCC 815 




12:25:37.1 


+ 13:08:36 


dE2,N 


16.1 


-700+50 


1 


VCC 817 


IC 3313 


12:25:37.7 


+ 15:50:06 


dEl 


15.0 


1086+57 


3 


VCC 823 




12:25:38.8 


+ 12:18:49 


dE0,N 


15.7 


1691+33 


2 


VCC 833 




12:25:44.4 


+ 13:01:20 


dE0,N 


17.48 


720+50 


1 


VCC 846 




12:25:50.3 


+ 13:11:52 


dEl,N 


16.2 


-730+50 


1 


VCC 854 




12:25:55.6 


+ 12:46:11 


dE8,N 


17.3 


684+75 


3 


VCC 856 


IC 3328 


12:25:57.9 


+ 10:03:14 


dEl,N 


14.42 


972+32 


2 


VCC 871 




12:26:05.6 


+ 12:33:34 


dE,N 


15.4 


1636+76 


3 


VCC 872 




12:26:00.7 


+ 12:51:40 


dEO.N 


17.0 


1183+50 


3 


VCC 882 


NGC 4406B 


12:26:15.6 


+ 12:57:51 


dE3,N 


16.7 


1101+55 


10 


VCC 896 




12:26:22.5 


+ 12:47:00 


dE3,N 


17.8 


2114+63 


3 


VCC 916 




12:26:33.3 


+ 12:44:24 


dEl,N 


16.2 


1115+32 


3 


VCC 917 


IC 3344 


12:26:32.4 


+13:34:43 


dE6 


14.8 


1375+66 


2 


VCC 928 




12:26:39.6 


+ 12:30:48 


dE6,N 


16.2 


-254+54 


3 


VCC 940 


IC 3349 


12:26:47.1 


+ 12:27:15 


dEl,N 


14.78 


1563+57 


9 


VCC 951 


IC 3358 


12:26:54.4 


+ 11:40:06 


dE2p,N 


14.23 


2066+21 


2 


VCC 953 




12:26:54.6 


+ 13:33:57 


dE5?p,N 


15.70 


-629+65 


2 


VCC 965 


IC 3363 


12:27:02.9 


+12:33:37 


dE7,N 


15.40 


790+50 


9 


VCC 967 




12:27:04.0 


+ 12:51:54 


dE4 


18.7 


1135+51 


3 


VCC 977 




12:27:11.1 


+ 12:02:18 


dE4,N 


17.9 


102+45 


3 


VCC 990 


IC 3369 


12:27:17.1 


+ 16:01:30 


dE4,N 


14.81 


1727+34 


2 


VCC 996 




12:27:21.2 


+ 13:06:36 


dE5 


18.4 


-28±32 


3 



Morphological types are from the estimates given by Binggeli, Sandagc & Tammann (1985) 
B (l) Binggeli, Popescu & Tammann (1993), (2) Binggeli, Sandagc, & Tammann (1985), (3) This Paper, (4) Young & Currie (1995), (5) Grogin et al. (1998), (6) 
Duprie & Schneider (1996), (7) Zabludoff & Mulchacy (1998), (8) Drinkwater ct al. (1996), (9) do Vaucouleurs et al. (1993), (10) Strauss et al. (1992), (11) Huchra 
ot al. (1993), (12) Rubin ct al. (1999), (13) Hoffman et al. (1996) 
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TABLE 3 

Virgo Dwarf Elliptical Galaxies with Radial Velocities Cont. 



Name 


ID 


R..A.(J2000.0) 


Dec. (,12000.0) 


Type* 


Magnitude 


vfkm/ s) 


Source^ 


VCC 1027 




12:27:38.0 


+ 12:52:48 


dE0,N 


18.1 


77±46 


3 


VCC 1036 


NGC 4436 


12:27:41.6 


+ 12:18:59 


dE6,N 


14.03 


1163+50 


2 


VCC 1040 




12:27:44.5 


+12:58:55 


dE3,N 


17.5 


-143+63 


3 


VCC 1042 




12:27:45.8 


+ 12:52:19 


dE4 


18.6 


-62+13 


3 


VCC 1073 


IC 794 


12:28:08.6 


+ 12:05:36 


dE3,N 


14.23 


1899+19 


2 


VCC 1075 


IC 3383 


12:28:10.7 


+ 10:17:43 


dE4,N 


14.95 


1844+40 


2 


VCC 1087 


IC 3381 


12:28:15.1 


+ 11:47:23 


dE3,N 


14.42 


645+27 


2 


VCC 1104 


IC 3388 


12:28:27.9 


+ 12:49:24 


dE5,N 


15.31 


1704+31 


2 


VCC 1122 


IC 3393 


12:28:41.7 


+12:54:57 


dE7,N 


14.82 


436+29 


2 


VCC 1164 




12:29:07.8 


+ 09:26:31 


dE6,N 


16.8 


1040+10 


9 


VCC 1173 




12:29:14.7 


+ 12:58:42 


dE5,N 


16.06 


2468+68 


8 


VCC 1185 


M87DW07 


12:29:23.4 


+ 12:27:02 


dEl,N 


15.56 


500+50 


1 


VCC 1254 


NGC4472DW08 


12:30:05.3 


+08:04:29 


dE0,N 


15.0 


1350+29 


2 


VCC 1261 


NGC 4482 


12:30:10.4 


+ 10:46:46 


dE5,N 


13.66 


1850+30 


2 


VCC 1308 


IC 3437 


12:30:45.8 


+ 11:20:34 


dE6,N 


15.1 


1721+45 


2 


VCC 1348 


IC 3443 


12:31:15.7 


+ 12:19:54 


dE0p,N 


15.64 


1679+39 


1 


VCC 1351 




12:31:17.5 


+13:49:42 


dE4 


16.0 


187+48 


3 


VCC 1355 


*IC 3442 


12:31:20.0 


+ 14:06:53 


dE2,N 


14.31 


1332+63 


11 


VCC 1386 


IC 3457 


12:31:51.3 


+ 12:39:21 


dE3,N 


14.43 


1426+60 


1 


VCC 1407 


IC 3461 


12:32:02.7 


+ 11:53:25 


dE2,N 


14.82 


919+55 


1 


VCC 1420 


IC 3465 


12:32:12.2 


+ 12:03:43 


dE4,N 


16.2 


1022+87 


2 


VCC 1431 


IC 3470 


12:32:32.4 


+ 11:15:46 


dE0,N 


14.51 


2025+75 


2 


VCC 1453 


IC 3478 


12:32:44.2 


+ 14:11:46 


dE2,N 


14.34 


1949+26 


2 


VCC 1489 


IC 3490 


12:33:13.8 


+ 10:55:44 


dE5,N 


15.89 


80+50 


2 


VCC 1491 


*IC 3486 


12:33:14.0 


+ 12:51:28 


dE2,N 


14.8 


903+42 


2 


VCC 1514 




12:33:37.7 


+ 07:52:16 


dE7,N 


15.1 


538+72 


2 


VCC 1539 




12:34:06.8 


+ 12:44:30 


dE0,N 


15.4 


1390+50 


1 


VCC 1549 


IC 3510 


12:34:14.7 


+ 11:04:18 


dE3,N 


14.63 


1357+37 


1 


VLU loot 


IC 3518 




+oy.o / .za 






1 i i n 1 K K 




VCC 1661 




12:36:24.0 


+ 10:23:00 


dE0,N 


15.7 


1400+50 


1 


VCC 1669 




12:36:30.4 


+ 13:38:16 


dE6,N 


16.2 


600+50 


1 


VCC 1713 




12:37:29.0 


+04:45:02 


dE 


15.1 


1655+25 


5 


VCC 1743 


IC 3602 


12:38:06.6 


+ 10:05:01 


dE6 


15.1 


1279+51 


1 


VCC 1826 


IC 3633 


12:40:11.2 


+09:53:46 


dE2,N 


14.87 


2033+30 


2 


VCC 1828 


IC 3635 


12:40:13.4 


+ 12:52:29 


dE2,N 


14.9 


1517+57 


2 


VCC 1857 


IC 3647 


12:40:53.2 


+ 10:28:33 


dE4,N? 


14.33 


634+69 


2 


VCC 1861 


IC 3652 


12:40:58.4 


+ 11:11:03 


dE0,N 


14.37 


683+54 


2 


VCC 1870 




12:41:15.3 


+ 11:17:55 


dE6 


15.8 


1617+27 


3 


VCC 1876 


IC 3658 


12:41:20.4 


+ 14:42:02 


dE5,N 


14.85 


45+49 


2 


VCC 1881 




12:41:29.5 


+ 10:45:46 


dE3;,N 


16.2 


138+56 


3 


VCC 1886 




12:41:39.4 


+ 12:14:52 


dE5,N 


14.87 


1159+65 


2 


VCC 1890 


IC 3665 


12:41:46.3 


+ 11:29:15 


dE4p 


14.84 


1227+57 


2 


VCC 1891 




12:41:48.9 


+ 11:11:29 


dE4,N 


16.7 


1016+43 


3 


VCC 1897 


UGC 7857 


12:41:54.2 


+ 13:46:22 


dE4,N 


14.49 


18+60 


2 


VCC 1910 


IC 809 


12:42:07.8 


+ 11:45:16 


dEl.N 


14.17 


206+26 


2 


VCC 1919 




12:42:19.0 


+ 10:34:04 


dE0,N 


17.0 


1869+54 


3 


VCC 1941 




12:42:50.0 


+ 13:17:38 


dE 


18.0 


1213 


12 


VCC 1947 




12:42:56.3 


+03:40:35 


dE2,N 


14.56 


944+35 


5 


VCC 1948 




12:42:58.0 


+ 10:40:55 


dE3 


15.1 


1581+22 


3 


VCC 1949 


NGC 4640 


12:42:58.1 


+ 12:17:10 


dE6,N 


14.19 


2077+75 


9 


VCC 1958 




12:43:10.0 


+ 11:02:12 


dE2,N 


17.0 


1049+30 


3 


VCC 1971 




12:43:31.0 


+ 11:02:50 


dE3: 


16.6 


1376+34 


3 


VCC 2012 


IC 3727 


12:45:05.7 


+ 10:54:03 


dEl:,N 


14.30 


2230+11 


9 


VCC 2019 


IC 3735 


12:45:20.4 


+ 13:41:33 


dE4,N 


14.55 


1895+44 


2 


VCC 2050 


IC 3779 


12:47:20.7 


+ 12:09:59 


dE5,N 


15.2 


1367+75 


2 


VCC 2062 




12:47:59.9 


+ 10:58:33 


dE 


19.0 


1146+8 


13 



Morphological types are from the estimates given by Binggcli, Sandagc &c Tammann (1985) 
B (l) Binggcli, Popescu & Tammann (1993), (2) Binggeli, Sandagc, & Tammann (1985), (3) This Paper, (4) Young & Currie (1995), (5) Grogin et al. (1998), (6) 
Duprie & Schneider (1996), (7) Zabludoff & Mulchaey (1998), (8) Drinkwater ct al. (1996), (9) do Vaucouleurs ct al. (1993), (10) Strauss et al. (1992), (11) Huchra 
ct al. (1993), (12) Rubin et al. (1999), (13) Hoffman et al. (1996) 
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TABLE 4 

Mean and g for Virgo "Core" Populations 



Type 


Mean (km s v ) 


a (km s x ) 


Median (km s x ) 


N 


E 


1090±73 


462 


1095 


40 


SO 


1164±86 


647 


1115 


56 


dE 


1146±150 


636 


1220 


18 


dE,N 


1033±84 


747 


1159 


79 


dET 


1054±73 


726 


1159 


97 


dSO 


960±116 


612 


895 


28 


Irr 


1054±81 


727 


1100 


81 


S 


1056±71 


776 


1037 


119 


Virgo 


1064±34 


705 


1100 


429 
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TABLE 5 

K-S Test Results for Radial Velocities 





E 


SO 


dE 


dE,N 


dET 


dSO 


Irr 


Spirals 


E 


100 


38.7 


60.9 


39.1 


36.8 


29.7 


36.8 


9.6 


SO 


38.7 


100 


62.0 


34.7 


48.0 


26.7 


59.1 


19.7 


dE 


60.9 


62.0 


100 


66.1 




46.0 


53.5 


35.4 


dE,N 


39.1 


34.7 


66.1 


100 




70.0 


85.4 


67.3 


dET 


36.8 


48.0 
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TABLE 6 

K-S Test Results for Projected Radial Positions 
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TABLE 7 

Projected Spatial and Velocity Gradient Slopes 



Type 



dlnJ dlno-^ 
dlnR dlnR 



Ellipticals -2.04±0.85 -0.45±0.10 

Irregulars -1.01±0.22 -0.08±0.11 

dETs -1.52±0.10 -0.09±0.05 

Spirals -2.01±0.30 -0.06±0.10 



-34- 




4000 5000 6000 7000 



Wavelength (A) 



Fig. 1. — A typical WIYN Hydra Multi-Object Spectrograph spectrum of a dwarf elliptical. 
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Fig. 2. — The difference between previously derived radial velocities and values measured in this study. 
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Fig. 3. — The projected distribution (J2000) of Virgo galaxies with known radial velocities. The large X 
symbol represents the location of the large elliptical M87. The 6° circle represents the spatial area used to 
define the Virgo core sample. The different points for galaxies represent different velocities. The larger the 
symbol, the larger the difference from 1000 km s _1 , with triangles for V< 1000 km s _1 , and circles for 
V> 1000 km s- 1 . 
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Fig. 4. — Velocity histograms of six galaxy populations contained in the entire Virgo sample. The solid part 
of the dE and dE,N histogram represents the non-nucleated dEs, while the shaded area is the combined dE 
and dE,Ns. 
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Fig. 5. — Velocity histograms for the Virgo 'Core' population, with fitted Gaussian profiles. The galaxies 
here are the subsample of galaxies in Figure 4 that are within 6° from the dynamical center of Virgo (sec 
text) and whose radial velocities are < 2400 km s _1 . The solid part of the dET histogram represents the 
pure dEs. 
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Fig. 6. — Positional plots of the various Virgo populations. These contain all galaxies considered by Binggeli 
ct al. (1985) to be within Virgo based on velocity and/or morphological information. The X symbol 
represents the location of M87. 
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Fig. 7. — Relative velocities of Virgo galaxies plotted as a function of their absolute magnitudes assuming 
a Virgo distance modulus of (m-M) = 31.3. From this plot there is no sign of mass- velocity segregation or 
energy equipartition in Virgo. The dashed lines represent the escape velocity of Virgo. 
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Fig. 8. — Elliptical, spiral and dET velocities as a function of magnitude. The open stars represent the 
velocity dispersion as a function of magnitude. The two upper X symbols in the elliptical galaxy plot arc 
two galaxies in our sample that are likely not part of the Virgo cluster proper or true ellipticals (see text). 
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Fig. 9. — Declination and right ascension velocity slices through the core of Virgo showing the 'finger of god' 
signature and infall ring patterns. The spirals arc X symbols, ellipticals black squares and dwarf ellipticals 
open circles. Several other clumpy regions can be seen in this diagram, most particularly at (-3° , 1000 
km s- 1 ) and (-5° , 2000 km s" 1 ). 
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Fig. 10. — The radial velocities of the Virgo galaxy populations plotted as function of projected distance 
from the center of Virgo. The high velocity dispersions (open stars) at the center of Virgo indicate that 
galaxies in Virgo are on radial or highly elliptical orbits. The black squares represent the average velocity 
at each distance. 
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Fig. 11. The ratio of the increase in internal energy produced in an impulsive interaction to the total 
internal energy, for two model galaxies as a function of the impact parameter b. For both a galaxy with a 
large size (5 kpc: solid line) and a smaller size (2 kpc: dashed line) the increase in internal energy is not 
enough to destroy the galaxy with impact parameters greater than <~ 30 kpc. 
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Fig. 12. — The projected spatial distribution, N, in galaxies per square degree as a function of radius. 
Ellipticals have the lowest densities and are represented by a solid line. The spirals are short-dashed lines 
and irregulars long-dashed ones. The dETs, the densest population, are designated by the small dotted line. 



